Experimental disagreements with Extended Unitarity 
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Abstract 

In production processes, e.g. J/^' wvrvr or pp — > Sir, the a and /o(980) overlap in 

the same partial wave. The conjecture of Extended Unitarity (EU) states that the tttt 

pair should have the same phase variation as vrvr elastic scattering. This is an extension of 

Watson's theorem beyond its original derivation, which stated only that the s-dependence 

of a single resonance should be universal. The prediction of EU is that the deep dip 

observed in vrvr elastic scattering close to 1 GeV should also appear in production data. 

^ Four sets of data disagree with this prediction. All require different relative magnitudes of 

O ■ a and /o(980). That being so, a fresh conjecture is to rewrite the 2-body unitarity relation 

, for production in terms of observed magnitudes. This leads to a prediction different to 

' EU. Central production data from the AFS experiment fit naturally to this hypothesis. 

^ ! 

. PACS numbers: 13.25.-k, 13.25. Gv, 13.75.Lb 



1 Introduction 



g_) ' In its simplest form, the idea of Extended Unitarity (EU) states that the vrvr pair in a single 
p ^. partial wave should have the same phase variation with s in all reactions as in elastic scattering. 
1^ ! This idea originates from Aitchison [1] and has been adopted in various guises by many authors. 
His arguments will be presented in detail in Section 2, so as to expose the assumptions and 
consequences. At the time the idea was introduced, it was a reasonable conjecture; now modern 
> ■ data allow it to be checked accurately, but disagree with it. 

^ ■ Many experimental groups have made extensive fits to production data using a /^-matrix 

Q\ ', approach based on EU. These fits are excellent; no criticism is intended of their quality. Exper- 
imentalists have found empirically the necessary freedom to get good fits to data. However, on 
close inspection, this freedom is inconsistent with strict EU. 

This whole topic has been the subject of extensive discussion with many authors. There is 
O ■ a bewildering jungle of claims and counter-claims. My objective is to cut a path through this 
tangle and expose where problems lie; this makes the presentation pedantic in places. 
■ Aitchison's essential point is that all processes should be described by a universal denominator 

?H I [1 — ip{s)K{s)], where K is the same as for elastic scattering; p is Lorentz invariant phase space. 
The assumption which is being made is that Watson's theorem [2] applies to the coherent sum of 
all components in the = 0^ partial wave. This is a step beyond Watson's derivation, which 
referred only to a single eigenstate; Watson did not consider overlapping resonances. 

In TTTT elastic scattering, the /o(980) is superimposed on a slowly rising amplitude associated 
with the a pole. Cern-Munich data [5] show that the phases of these two components add. 
Below the KK threshold, both a and /o(980) T-matrices T = e*^sin5 are confined to the 
unitarity circle if we neglect the tiny inelasticity due to vrvr — > 77. Unitarity may be satisfied 
by multiplying S'-matrices S = e^*^ of a and /o(980), as suggested by Dalitz and Tuan [3]. This 
fits the data within errors of ~ 3.5°. 

Fig. 1 shows the Argand diagram for the / = tttt S-wave from my recent re-analysis of 
these (and other) data [S]. From BES data on J/\l/ 07r+7r~, the /o(980) has a full-width 
at half-maximum of 34 ± 8 MeV [B]. The combined phase shift rises rapidly from 90° at 0.88 
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Figure 1: Argand plot of the tttt / = S-wave in elastic scattering; masses are marked in GeV. 



GeV to 270° near 1.1 GeV. There is a deep dip in the cross section where the combined phase 
goes rapidly through 180°. The crucial point of EU is that this feature should be common to 
production processes. 
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Figure 2: The tttt mass projection for BES data on J/\E' ^ ^tt+tt : the upper histogram shows 
the current fit to experimental points [5]; the lower histogram shows the fitted a component. 



BES data on J/\E' — ^vrvr [S] immediately require a modification of the rudimentary form of 
EU. The TTTT mass spectrum in these data is reproduced in Fig. 2. There is a dominant /o(980) 
contribution and a small interfering a contribution; this is very different to elastic scattering. 
Lahde and Meissner [7] modify the conjecture of EU to apply separately to strange and non- 
strange components, i.e. to the scalar form factors for tttt and KK. 

The dip in the elastic cross section at 989 MeV is a very delicate feature. If, for any reason, 
relative magnitudes of a and /o change, the zero at 989 MeV can disappear quickly; here and 
elsewhere, /o will denote /o(980) unless there is confusion with other /o's. If the phase of /o 
changes with respect to the a, the mass at which the dip appears will likewise change. The 
interference region between a and /o(980) is an ideal place to check Extended Unitarity. 

Two considerations will play a critical role: unitarity and analyticity. Unitarity is often 
quoted and plays an essential role in setting up the current ii"-matrix formalism which treats 
both elastic scattering and production on the same basis. For a production reaction, Aitchison 
conjectures a unitarity relation for the production amplitude F: 



lmF = FT*, 



(1) 



He defines his F to be proportional to T^^^ / p, where T^^^ refers to production: 



Fp 



Dividing both sides of (1) by a, 



IniT^f) 




(3) 



It is odd that T^^^ on the right-hand side is multiphed by T*i, unless T^^ = T^i. However, 
experiment will require different contributions to T^^^ from a and /q. 

Consider next analyticity. Dispersion relations connect magnitudes and phases. If the relative 
magnitudes of /o and a change from those of elastic scattering because of matrix elements, their 
relative phases must also change. Conversely, analyticity predicts that if the phase variation 
with s of the amplitude is universal, as EU demands, so is the variation with s of the magnitude 
(up to a constant scaling factor); for the simplest situation where only resonances are present, 
the relative magnitudes of a, /o(980) and any further /q must be almost the same in production 
as elastic scattering. This is a point which has almost always been ignored. 

The word 'almost' represents a caveat: there may in addition be a polynomial is s which 
can be different between elastic scattering and production. It turns out that one can plausibly 
limit deviations of relative magnitudes within 12%. This question is discussed in subsection 
2.1. Experiment requires larger deviations than this in the four sets of data discussed here. 
This implies phases must change from those predicted by EU. Experimentalists have correctly 
allowed for this by using complex coupling constants in the isobar model. 

Section 3 compares the prediction of EU with 3 sets of data. The first concerns BES data for 
J/\I' — > tUTT+TT^ [8j. The non-strange components of cr and /o dominate both elastic scattering 
and J/\E' utttt. From Aitchison's algebra and that of Lahde and Meissner, it follows that 
the /o amplitude should have almost the same magnitude as the a amplitude, as well as the 
same phase as elastic scattering. This prediction is contradicted by the data, where no /o(980) 
is visible and a fit to the data places a low limit on it. 

The next two sets are Crystal Barrel data for pp 37r°, where a and /o(980) are clearly 
visible [9]. One set is for annihilation in liquid hydrogen and the other for gaseous hydrogen. 
Annihilation from the ^Pi initial state is 13% in liquid and 48% in gas, allowing a clear separation 
of amplitudes for production of a and /o from ^5*0 and ^Pi. Results for both are inconsistent 
with the deep dip of elastic scattering predicted by EU. 

Section 4 concerns data from the AFS experiment on central production: pp ppvT+vr" [lOj . 
Here one expects the protons in the final state to act as spectators. However, EU still fails 
conspicuously to fit the data. This important result leads to a revised form of the unitarity 
relation, as follows. 

Fig. 3 sketches the usual diagrammatic approach to the unitarity relation. It may be derived 
by cutting the diagram down the middle, along the dashed line. For a 2-body system of tttt, 
KK, riT], etc. the resulting relation is well known: 



The application of 2-body unitarity assumes that the pions interact only with one another, not 
with any spectator. In most sets of data there are large signals where pions do interact with the 
spectator. For J/\Ef — *• ujtttt, as an example, the 6i(1235)7r channel accounts for 40% of events 
[H]. Some of it may be generated by pions from decays of a or /q rescattering from the spectator; 
this is a so-called triangle graph. Aitchison himself remarks that this can distort the unitarity 
relation. This provides one reason why EU may fail for the first three sets of data; it does not 
explain the fourth, where some further effect is required. 



ImTf 



el - J^eUel- 



(4) 
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Figure 3: Unitarity diagram for tttt — > tttt. 



There are fundamental differences between elastic scattering and production. In J/\E' —>■ 
UTTTT, for example, matrix elements < J/'^lua > and J/\E'|co'/o > dictate the magnitudes of 
these amplitudes; any values are possible. This differs from elastic scattering, where a and /o 
magnitudes are fixed purely by their coupling constants g.,, to vrvr. Equation (3) is an asymmetric 
relation, allowing a and /o to be produced with different magnitudes, but requiring that they 
rescatter an in elastic scattering. A more logical alternative is the symmetric relation 

I^T'^p) = T^p)t(p> , (5) 

hence ImF = FT*^'p\ This relation fits AFS data for central production naturally, whereas 
EU does not. If the /o is absent from production data, (5) reduces to the obvious relation 
ImT*^ = iT'^p; Aitchison's form of the relation, taken with analyticity does not allow the /o to 
be absent, as we shall see in Section 2. 

Section 5 suggests a new way of fitting 2-body data. Section 6 then summarises conclusions. 



2 The hypothesis of Extended Unitarity 

In a two-body process, the scattering of a pair of pions to final states nn, KK, rjrj, Att and 77 
must satisfy unitarity. The T-matrix for these coupled channels may be written in terms of a 
real i^T-matrix as 

T,i = pKil-tpK)-' . (6) 
It is normalised here so that T^tt = (776^**^ — l)/2i. Below the inelastic threshold 

pK = tan 6. (7) 

The T-matrices used here will include couplings to all channels. However, it simplifies the 
presentation of essential points to reduce the formalism initially to a single tttt channel. This 
simplication is sufficient to expose the basic issues, and can be generalised later to include 
inelasticity. 

The approach of Aitchison [T] will now be outlined. I am grateful to him for clarifying the 
algebra in more detail than is to be found in the original publication. Suppose S-matrices 



multiply, i.e. phases add. Let Ka and Kb be i^-matrices for a and /o respectively. The 
elementary expression for tan((5A + Sb) then gives a X-matrix for elastic scattering 



- 1-p^KaKb^ 



from which it follows that the T-matrix for elastic scattering is 

{l-ipKA){l-ipKB)' 



'el 



(9) 



Aitchison now conjectures that an amplitude F for producing a two-body channel present in 
K^i may be written in terms of a vector P, with 

F = {l-ipK,i)-'P (10) 

_ aAKA + OiBKB ,^^s 

^ " I-P^aKb ' ^ ^ 

where a a and are constants for production couplings. With this ansatz, the relation 

ImF^FT*, (12) 

known as Extended Unitarity, is automatically satisfied. It is a consequence of (10) that F has 
the same phase as T^i. Substituting (8) and (11) in (10) gives, in this one-channel case 

{1-ipKa){1-ipKb) ^ ' 

= aA[TA{l + iTb) + I3Tb{1 + iTa)]/ p, (14) 

where = aB/o-A- From (13), the phase of F is indeed 5a + ^s, as imposed by (10). Equation 
(14) will play the decisive role in comparisons with experiment. 

In (14), TAi^+iTs) = expi(5^+5s) sin^y^cos^B andTB(l+iT4) = cxpi{6A+SB) cos ^4 sin 5^. 
At 989 MeV, 6f = 90° and = 92°. So both terms are very close to zero, regardless of the 
values of a a and a^. This predicts that production data should have the same deep dip at this 
energy as elastic scattering. 

There is a further point. In the second term, (1 + iTa) — over a sizable mass range. In the 
first term, Tb should be conspicuous, since it has a rapid phase variation and the same peak 
magnitude as Ta, which is itself clearly visible in all sets of data considered here. However, the 
data all require the magnitude of the /o(980) to be smaller than predicted. 

The key point is that the factor 1/(1 — ipKb) of (13) leads directly to the factor (1 -|- iTb) in 
the first term of (14). The first and third sets of data will require fi of (14) to be small. In the 
elastic region, the first term becomes 

F ^ iaA{l + ie'^" Sin5B). 

The bizarre conclusion of EU is that the phase of the /o(980) is present even though aB — 0. 
This is inconsistent with analyticity. It will be shown in Section 3 that experiment disagrees 
with EU even without the constraint of analyticity. However this additional constraint makes 
conclusions more definitive. 
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2.1 Analyticity 

For purely elastic scattering, the Omnes relation [11] reads [including a factor p{s) in A^(s)]: 



N{s)/D{s) 



(15) 



e 



i5{s) 



exp — 




' — Ami s — s' 



ds' 5{s') 



(16) 



where P denotes the Principal Value integral. We shall not actually need to evaluate (16). It 
plays only a conceptual role and this needs considerable explanation. The basic point is that 
D[s) contains both real and imaginary parts, so 5{s) determines both. For elastic scattering 
N{s) is real. It arises from the left-hand cut, i.e. meson exchanges between the two pions. 

With inelasticity, corresponding relations may be written in a 2-channel form. Then 5{s) 
is replaced by 0(s), the angle T makes to the real axis when measured from the origin of the 
Argand diagram, see Fig. 12 of subsection 4.1. 

If EU is valid, the production amplitude may be written F = X{s)/D{s). In principle 
X{s) could be anything, depending on production dynamics. However, we have quite precise 
experimental information about it. An extreme view is that and as of (13) can be arbitrary 
and complex. However if they are complex this leads directly to a conflict with EU. Equation 
(14) contains two parts Ta{1 + iTb) and Tb(1 + iTa)- Substituting Breit-Wigner formulae for 
Ta and Te, the first term becomes g\{M'^ — s) / Da{s)Db{s). This is real but has a specific 
s-dependence in the numerator as well as in the denominator. If a a or becomes complex, 
the numerator becomes complex. This then introduces a phase variation separate from D{s). 
The "prediction" of EU is distorted by this extra phase. Only if X{s) is real does EU survive 
in its strict form. 

Many experimental groups have used the P-vector approach using complex coupling coeffi- 
cients, without realising that this destroys the universality of the phase variation with s. This 
is what experiment demands, so they have done the right thing. But the use of the universal 
denominator [1 — ipK{eiy\ is no longer logically correct. One might as well fit directly in terms 
of complex coupling constants and individual T-matrices for each resonance. 

2.2 Form Factors 

There is a further fundamental point. For elastic scattering, N{el) is uniquely related to ImD(s) 
by both unitarity and analyticity. At first sight it appears that analyticity relates X{s) in the 
same way to D{s) in production reactions, with the result X{s) = aN{el), where a is a constant. 
This requires /? = 1: if the phase of the tttt amplitude is universal, relative magnitudes of a and 
/o must also be universal. 

There is however a caveat. A more fundamental form of the unitarity relation (12) is that 
the discontinuity of F across the elastic branch cut is 2iFT*. Then F may be multiplied by 
a polynomial X{s), providing it does not have a discontinuity along the real s-axis. A few 
examples will hopefully clarify ideas. Firstly, a form factor in s is one such example, arising 
from the sizes of particles, i.e. from matrix elements. Secondly, in — »• 7/0, the El transition 
has an intensity proportional to the cube of the photon momentum; this infiates the lower side 
of the /o(980). Thirdly, in ^Pi pp — > tt^ct, there is an L = 1 centrifugal barrier for the production 
process. Fourthly, in some special cases, matrix elements may go through zero as a function of 
s. Taking X{s) to be real, let us write in general 



F = X{s)N{el)/D{s). 



(17) 
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A feature of all production data considered here is a strong low-mass vrvr peak due to the 
a pole. This peak is not present in elastic scattering because of an Adler zero in the elastic 
scattering amplitude at sa — 0.41m^, just below the tttt threshold. The elastic amplitude rises 
approximately linearly with s and there is no low mass peak. The origin of the difference has 
been known to theorists for at least 20 years. Au, Morgan and Pennington [12] pointed out 
that the difference between elastic scattering and central production data can be accomodated 
by using the same Breit-Wigner denominator for both, but replacing the numerator N{el) by 
something close to a constant. This polynomial is allowed because sa is outside the physical 
region. Data require X{s)N{el) ~ 1, hence X{s) ~ 1/N{el) ~ — sa)- More exactly, 

X{s) = l/[{s - sa){1 + bs) exp[-(s - M^/A] (18) 

for the parametrisation of the cr amplitude in [13]. In practice, quadatic and cubic terms in s 
are very small and under tight control from fits to data up to 1.8 GeV. 

For J ^ uTiTi, the a pole is visible by eye in Fig. 4(a) below. The phase of the a amplitude 
in this reaction is experimentally the same as in elastic scattering within ~ 3.5° [14J. Values 
of N{prodn) = X{s)N{el) can be determined directly from the data. The same is true of the 
K [15], which likewise has an Adler zero in the numerator for elastic scattering, but not for 
production. In both cases, N{prodn) is consistent within errors with a constant; the Adler zero 
in the numerator of elastic scattering has disappeared. One can try fitting the a and k poles 
in production data with the conventional form factor N{prodn) = exp(— fc^i?^/6), where k is 
momentum in the production channel. For both, optimises at slightly negative values, which 
are unphysical. For the a, < 0.30 fm^ with 95% confidence and for the k, R^ < 0.38 fm^ at 
the same level. 

A crucial piece of information in testing EU will be relative magnitudes of a and /q. The 
magnitude of the a amplitude is easily separated from the tail of the /o(980) at 920 MeV, 
three half- widths from 989 MeV. If X{s) is determined at this energy and at 400 MeV, the a 
amplitude changes by 20% at 989 MeV for the 95% confidence level quoted above. However, 
this exaggerates the error, since there are compensating changes in the line-shape fitted to the 
amplitude. Realistically, changes are half this. Adding in quadrature uncertainties in the a line- 
shape due to the opening of the KK threshold, the uncertainty in the amplitude extrapolated 
from 920 to 989 MeV is < 12% with 95% confidence. This provides a tight constraint on the 
relative magnitudes of a and /o(980) amplitudes if EU is correct. This disposes of the paradox 
that the phase of /o(980) can be present with as = 0. 

From this point onwards, it will be assumed that EU should be supplemented with the 
condition imposed by analyticity within 12%. 

2.3 Formulae for a and /o 

Formulae and numerical parameters for the a amplitude are to be found in Refs. [13] and [5]. 
The TTTT coupling has been fitted to four sources: (i) phase shifts deduced from Cern-Munich 
data by Ochs [16j, (ii) Ke4 data of Pislak et al. [T7], (iii) predictions of nn phase shifts by 
Caprini et al. [l^ using the Roy equations, and (iv) BES data on J/\E' — > cutt+tt^ [8]. 

The coupling to KK and rjr] has been fitted to available data on those channels [19], and in 
[5] the 47r coupling has been derived from Cern-Munich data. There is close consistency between 
all these sets of data. Because Refs. [13] and [5] fit the same data with different formulae, the 
amplitudes agree within errors up to 1.2 GeV. Those of Ref. [S] are more cumbersome to use, 
since they allow for the dispersive effect of the Air threshold. Therefore the first and fourth sets 
of data discussed below are fitted with the formulae from Ref. [13] . 
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The general procedure adopted here is to allow the parameters of the a to have the freedom 
allowed in earlier determinations of its parameters, but no extra freedom in the vicinity of 
/o(980). In testing EU, the magnitude of /o is restricted to the 12% discussed above; its phase 
is allowed the freedom with which its parameters are known. 

2.4 Parameters of /o(980) 

The /o(980) is so narrow that it is readily separated from the a. The BES data on J/\l/ (j)7T~^7T~ 
shown in Fig. 2 exhibit a very clear /o(980) —>■ vrvr signal. An important point is the excellent 
mass resolution and mass calibration of the BES detector, ~ 4 MeV. Both are easily checked 
for the KK channel against the very precisely known parameters of 0(1020). Data from the 
same publication [6J on (f)K^K~ contain a clear /o(980) —>■ KK peak, and the two sets of data 
determine accurately the ratio g'^{KK)/ g'^{TnT) of couplings to KK and vrvr. An important detail 
is an error in units in [6]: g'^{7nr) is reported as 165 MeV; this should read 0.165 GeV^. 

The Breit-Wigner denominator for the /o(980) amplitude is [M^ — s — i{glpi + ^72^2)] and 
P2 has to be continued analytically below the KK threshold as i^jAmj^/s — 1. Without direct 
information on the KK channel, this term gets confused with (M^ — s) [20]. Any form factor 
applied to g'^{KK) adds to the confusion. One only has to glance at the Particle Data Tables [21] 
to see the large spread in parameters fitted to the /o(980) (and ao(980)) in experiments having 
no direct information on the branching ratio between KK and vrvr. Unfortunately, the PDG 
does not report the BES determination of g'^{KK)/g'^{TnT), which is the best in the published 
literature because of the availability of a clear signal in (pKK. 

It was shown in [19] that the BES parameters for /o(980) are closely consistent with Kloe data 
[22j on — > vr°vr°7 when one allows for interference between a and /q. This paper determines 
g'^{a KK)/g^{a vrvr) = 0.6±0.1. BES data on J/^ ujK+K~ confirm a large value > 0.6 
for this ratio [23j . The vrvr full width at half- maximum (34 MeV) agrees well with Cern- Munich 
data 30 ± 10 MeV [3]. It also agrees closely with the full-width of the /o(980) signal in Crystal 
Barrel data (~ 46 MeV). The BES parametrisation will therefore be adopted in fitting the AFS 
data. Further checks from Belle, Babar and Cleo C will be very welcome. 

3 Experimental tests 

3.1 J/^' U7T7T 

Considerable detail needs to be given of fits to experimental data, in order to pin down the 
disagreements with EU. The prime conclusions which will emerge are that (i) the vrvr amplitude 
does not follow that of elastic scattering, (ii) the magnitude of the /o(980) amplitude, relative 
to a, is much smaller than predicted by Eq. (14), regardless of analyticity which also requires 
that their relative magnitudes should be the same within 12%. 

Fig. 4 displays features of BES data for J/^ cuvrvr [8]. There is a clear peak in Fig. 4(a) at 
~ 0.5 GeV due to the a pole. Its shape is cleanly separated from the /2(1270) contribution up to 

1.05 GeV, where the /2 rapidly overtakes it. The data include a slowly varying 14% background 
which is included in the fit. There is also a well defined slowly varying component due to the 
reflection of 61 (1235). Both are shown in the experimental publication. 

There are two amplitudes for production of a and /o(980), with orbital angular momenta 
L = and 2 in the production reaction. The L = 2 amplitude includes a centrifugal barrier for 
production which optimises at a radius of 0.8 fm (roughly as expected for convolution of wave 
functions of a, uj and /o). For the isobar model fit, different magnitudes are allowed for L = 
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Figure 4: Fits to BES data for J/ip ^ uj7t~^tt~: (a) vrvr mass projection from the full Dalitz plot 
using the cr amplitude of Ref. [13]; the lower histogram shows the S-wave component; (b) Dalitz 
plot from data; (c) as (a) using Eq. (14) for EU; (d) enlarged view of part of the Dalitz plot. 



and 2, with complex coupling constants, though it turns out that the L = 2 contribution from 
/o(980) is negligible. In the EU fit the relative phases of a and /o components are constrained 
to be the same, and relative magnitudes are constrained to be the same within 12%. In Figs. 
4(a) and (c), the lower histogram shows the = O"*" intensity: from the isobar model fit in (a) 
and from EU in (c). The upper histogram, close to data points, shows the result of the full fit. 
Unfortunately, the mass projection alone is not definitive, for reasons explained shortly. One 
low point in the mass projection just above 1 GeV hints at a dip following the EU prediction. 

However, there is a great deal more information contained in the Dalitz plot of Fig. 4(b) and 
also in the correlation between the u decay plane and the utinT production plane. In the Dalitz 
plot, there are strong horizontal and vertical bands due to 6i(1235). These bands interfere with 
a and /2(1270) (and /o(980) if present). Note that the 6i(1235) decays mostly through S'-waves, 
so its intensity would be almost constant across the Dalitz plot in the absence of interferences. 
Also note from the lower histogram of Fig. 4(a) that the a amplitude at 950 MeV is sizable. 
Cross-hairs on Fig. 4(b) show where a tttt pair of mass 1 GeV intersects the bi band. If /o(980) 
were present with the same magnitude as a and with the phase predicted by EU, a large but 
narrow interference with bi cannot be avoided due to /o(980), whatever the relative phase of 
6i(1235). It should have a full width of ~ 0.07 GeV^ along the bi band: the line-width of 
/o(980) —>■ TTTT. There should also be a dip somewhere along the diagonal at mTrvr = 1 GeV. Fig. 
4(d) shows an enlarged view of this region. There is no sign of these features. 

The LJ decay plays an important role. The spin of the lv lies along the normal to its decay 
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plane; information on this decay plane is a key ingredient in determining helicity amplitudes. 
These are essential to determine interferences between the amplitudes for lou^ U!f2{1270) and 
61(1235)77. There are in principle five /2(1270) amplitudes corresponding to production with 
L = (one amplitude), L = 2 (three) and L = 4 (one). It turns out that the L = 4 amplitude 
and one of the three L = 2 amplitudes are negligible. There are large interferences between a 
and the remaining three /2(1270) amplitudes and even larger interferences with 61 (1235). 

With the EU hypothesis, the best fit to the Dalitz plot (upper histogram) and the lo decay 
plane fills in the predicted dip at s^ttt = 1 GeV^ with other interferences, see Fig. 4(c); however, 
the price is a considerable deterioration of log likelihood compared with the fit of Fig. 4(a). 
The isobar model fit is better than EU by 259 in log likelihood. There are two reasons. Firstly, 
the narrow /o docs not appear in interferences with either 61 or /2. The mass resolution of 
the BES data is 4 MeV in tttt at 950 MeV; searching for the /o(980) with accurately known 
line-shape is limited purely by statistics and there are ~ 40K events. Secondly, in the fit based 
on EU, there are strong confiicts between interference terms amongst 6i(1235)7r, t(;/2(1270) and 
uia. The data want slowly varying interferences between bi and the broad /2, instead of rapidly 
varying interferences of the narrow /o with bi and/or /2. 

There are two extra fitting parameters for the isobar model fit: i.e. two complex coupling 
constant for uia instead of one. [The /o(980)a; L — 2 amplitude is negligible]. The definition 
of log likelihood is such that a change of 0.5 corresponds to a one standard deviation change in 
one variable. For large statistics, is approximately twice the change in log likelihood. So the 
fit of Fig. 4(a) is better than EU by 19.7 standard deviations statistically. 

It is important to remark that the BES publication gives a second fit done independently. 
This shows a tttt mass projection almost identical with the scalar form factor. This requires a 
very strong /o(980), in conflict with what one can see by eye in Fig. 4(d). However, this analysis 
made no use of the uj decay plane. Without that information, it is impossible to disentangle 
the magnitudes of the five /2(1270) amplitudes; the determination of the magnitude and phase 
of the /o(980) signal becomes very poor. I have rechecked the analysis omitting information 
from the uj decay plane. No stable fit emerges with the tttt mass projection of the second fit in 
the BES publication, and there is no evidence for the presence of /o(980) at all in this fit. The 
only explanation 1 can find of this second BES fit is that the relative magnitudes of /o(980) and 
a have been constrained to agree with the scalar form factor. This is, of course, exactly what 
needs to be checked in the present work. 

Fig. 5 shows the Argand diagram for T{prodn) — Fp from my fit to a and /o; the factor p is 
included so that the diagram goes to zero at the tttt threshold. The ratio of /o and a amplitudes 
is 0.12 ± 0.06 at 0.99 GeV and the /o lags the a by (42 ± 20)°. At the top of a large loop due 
to the a, there is a small structure due to /o(980). For comparison purposes, the dashed curve 
shows the a amplitude alone and the dotted curve shows the EU prediction normalised to the 
full curve at 0.47 GeV, the peak of the a. The essential point is that the deep dip in this latter 
curve near the /o mass is missing from the data. There is one qualitative feature of Fig. 5 which 
will be important later. The amplitudes in a production process arc free to add vectorially in 
any way required by the data. This situation is quite different from elastic scattering, where the 
amplitude is constrained to the unitary circle. 

It is now necessary to consider a number of systematic questions which might blur the argu- 
ment against EU. Corresponding remarks will apply to fits to the other three sets of data. 

Firstly, could the disappearance of the predicted /o(980) arise from effects of the KK thresh- 
old? The answer is no, because the /o appears clearly in elastic scattering and the prediction of 
EU is that it should be almost identical in production. A precise cancellation of the magnitude 
and phase of the predicted /o would be needed. This is ruled out by the known couplings of 
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Figure 5: Argand diagram for Fp; masses are marked in GeV. Full curve, fit; dashed curve, a 
alone; dotted curve, EU prediction 



both a and /o to KK. The analytic continuation for the /o(980) below the KK threshold is 
accurately determined over the small mass range concerned by BES data on ^vrvr and (pKK. In 
the a amplitude, the KK inelasticity rises over a mass scale of 200 MeV as shown in Fig. 11 
of Ref. [19]. Any flexibility in the analytic continuation of the KK term then has a scale of 
~ 300 MeV. Furthermore, it is closely constrained by the fit to Kloe data on — > 'jtt^tt^. So 
this explanation is highly implausible. 

Secondly, could the BES data be fitted assuming J/\l' ^ uKK, followed by KK tttt? 
This has been tested by adding T12 fitted freely. Its magnitude optimises at zero within errors. 
It improves log likelihood only by 2. 

3.2 An objection of Aitchison 

Thirdly, Aitchison argues [23] that a further term C might be added to the P-vector due to 
dynamics of the production process. A similar approach was used in fits by Bowler et al. and 
Basdevant and Berger to the ai(1260) to allow for the Deck effect ^25] [26] [27]. An additional 
term Ce^^^'^~^^^\cos 6a cos 6b — sin^yisin^s) appears in Fp. If C = iaA, the second term can 
cancel the term iTaTb of EU, leaving and the first of the additional terms, which is close to 
since 6a — 90°. This removes almost all the structure due to /o(980). However, this cancellation 
also leaves a term i exp + 6b) cos 6a cos 6b and when 6b is small this becomes i exp i6A cos 6a, 
which is much larger than the term Ta itself if C is a constant. This is ruled out by the data, 
so it becomes necessary to tailor the s-dependence of C to reproduce the magnitude of the tttt 
amplitude. 

However, this is not the end of the story. Phase information is also available. In [H] it is 
shown that interference between a and 61 (1235) measures the phase of the a amplitude (plus any 
background term C) and requires it to be the same as for elastic scattering down to 450 MeV 
with an error of 5° in the worst scenario, see Fig. 2 of that paper. It is worth mentioning that 
there is a similar result for the k in [TJ]. Fig. 4(a) of that paper shows that the k phase from BES 
data on J/^ ^ K^'k~K~'k^ agrees with the LASS effective range formula for elastic scattering 
down to 750 MeV within 3.5°. Fig. 4(f) of the same paper shows a corresponding agreement 
for E791 data on D — > Ktttt within a similar error. These results rule out any background 
different from a and k, with a magnitude larger than 12% and a phase difference above ±5°. 
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They of course agree with Watson's original statement that D{s) of a single resonance should 
be universal. 

There is no obvious source of the background proposed by Aitchison, as there was in [25], 
where a Deck background was visible in the data and made a 40% contribution to cancelling 
the ai(1260) amplitude. A similar fortuitous cancellation with unidentified backgrounds will be 
required also for all three of the following sets of data. In all cases, the additional background 
gives rise to a phase variation different to strict EU. If the background removes the /q, only the 
a is left, i.e. T^-, which contributes only part of the phase of elastic scattering; so the amplitude 
does not satisfy EU. 
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Figure 6: Argand diagrams for the nn S-wave of Ref. [5] for (a) ^5*0 annihilation, (b) '^Pi. 



3.3 pp 37T^ 

A fresh analysis of Crystal Barrel data on pp annihilation at rest to 3tt^ has been completed 
recently [5]. Data are available with statistics of ~ 700K events (and 0.5% experimental back- 
grounds) in both liquid and gaseous hydrogen; these allow a good separation of ^5*0 and ^Pi 
initial states. Although the /o(980) appears clearly in ^5*0 annihilation, the magnitude of the 
Tin S-wave amplitude does not go to zero on the Argand diagram near 1 GeV, as EU predicts. 
This Argand diagram is reproduced in Fig. 6(a). Its magnitude is smallest at 0.98 GeV, but is 
still very distinct from zero. For ^Pi annihilation, the fitted /o(980) amplitude, relative to a, is 
much smaller, see Fig. 6(b). The magnitude of the a amplitude is quite large near 1 GeV. 

The question arises how robust these solutions are. Could they 'bend' to accomodate EU? 
The published analysis requires some s-dependence of the numerator fitted to the a amplitudes, 
of the form (1 -|- As), with complex A. Can extra flexibility reach agreement with EU? A second 
point is that the analysis does not include the repulsive J = 2 vrvr S-wave. Could this bring 
conclusions into line with EU? 

The brief answer is definitely no, and will be presented graphically. Fig. 7(a) shows the vrvr 
mass projection for data in liquid from the current analysis; it fits the data points accurately. In 
Fig. 7(b), the a + /o combination of the isobar model is replaced by the EU combination (with 
the constraint from analyticity that magnitudes of a and /o(980) should be equal within 12%). 
Initially, only the coupling constant of this combination is refitted, leaving other amplitudes 
untouched; this is for the purpose of illustrating the change required by EU. The a component 
is large and cannot change much; a deep dip appears at 1 GeV because of the corresponding 
dip in the elastic amplitude. 

It is of course necessary to re-optimise all components. The resulting mass projection is shown 
in Fig. 7(c) and is still in severe disagreement with the data. A measure of the disagreement 
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Figure 7: tttt mass projection of data and fits: (a) from Ref. [5] in liquid hydrogen, (b) with the 
fitted EU amphtude but other amphtudes untouched, (c) after refitting all amplitudes, (d) as 
(c) in gaseous hydrogen, (e) and (f) with the 1 = 2 S-wave included in the fit. 



may be obtained from x^- Here, it is necessary to point out that even the fit of Fig. 7(a) has a 
larger than 1. This is probably because of the enormous statistics and small, slowly varying 
systematic errors in acceptance. The procedure adopted here is to renormalise to 1 per point 
for this fit and apply the same scaling factor to all other fits. The fit of Fig. 7(c) then has a 
renormalised of 40599 for 3500 bins; this is a 170 standard deviation discrepancy, allowing 
for the reduction in the number of fitted parameters by 2. 

A much better fit is possible if relative magnitudes of a and /o(980) are allowed to vary. The 
phase of the S-wave amplitude in Fig. 6(a) is close to EU, and only a 10 MeV shift is required 
for a perfect fit. This is consistent with the energy calibration and resolution of the Crystal 
Barrel detector. However, the fitted combination of amplitudes no longer agrees with the crucial 
equation (14). 

Fig. 7(f) shows the effect of including the 1 = 2 S-wave amplitude, using the formulae of 
Section 4.1 of Ref. [S]. There is only a rather small improvement, because the slow s-dependence 
of this amplitude cannot fill the narrow dip at 1 GeV. The renormalised f^^Hs to 34345, a 
discrepancy of ~ 130cr. 

It is not just the rnr mass projection which governs x^- One should inspect discrepancies in 
X^ all over the Dalitz plot. Fig. 8 makes such a comparison. Panel (a) shows discrepancies in 
X^ where the fit is above the liquid data and (b) the discrepancies where the fit is low. Panels 
(c) and (d) show results for data in gas. One sees striking systematic discrepancies all over the 
Dalitz plot, arising from interferences. Such discrepancies are almost completely absent from 
the fit of Fig. 7(a) using the isobar model. 

There are two further points for discussion. In the work shown on Figs. 7(e) and (f), the 
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Figure 8: Discrepancies in over the Dalitz plot for (a) and (b) liquid hydrogen, (c) and (d) 
gaseous hydrogen. In (a) and (c) the fit is above the data and in (b) and (d) below. 



1 = 2 S-wave was fitted without factoring out the Adler zero which occurs at s ~ — 0.41m^. 
If this is done (as for the a amplitude), the broad 1 = 2 amplitude gets confused with the a 
amplitude and leads to minor improvements all over the Dalitz plot. However, none of these is 
distinctive enough to require the 1 = 2 amplitude definitively. 

Secondly, could a more complicated polynomial than (1 + As) multiply the a amplitude and 
give a successful fit? Extensive tests were made in |5] with the objective of improving the 
fits reported there. If one chooses too free a polynomial, the fitted ^Pi component in gaseous 
hydrogen can fluctuate wildly from the 50% predicted from Stark mixing by Reifenrofer and 
Klempt [28]. To avoid this, the fitted ^Pi component is constrained within the range (50 ± 7)%. 
Unless the numerator of the EU amplitude is designed to include a narrow dip at 1 GeV, no 
large improvement is observed. 

One further observation from the Crystal Barrel data is worth reporting. Relative intensities 
of /2(1565) and /2(1270) are quite different in pp annihilation to those in elastic scattering. In 
pp data, the integrated intensities of these two resonances are equal within 12% after allowing 
for the (modest) effects of centrifugal barriers for production {L = 2 for ^Sq annihilation and 
L = 1 for ^Pi). However, in elastic scattering the fitted /2(1565) tttt width is a factor 4 
smaller than that of /2(1270) tttt. This again disagrees with EU. 



4 Central production of nn 

Central production of a tttt pair in pp pp{7T7i:) was fitted using EU by Au, Morgan and 
Pennington (AMP) [12]. In data from the AFS experiment at the ISR [10], the two final-state 
protons are produced with very small 4-momentum transfers t = —0.003 GeV^. It is routinely 
assumed that the tttt pair is unaffected by final state interactions with these protons, which are 
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separated from the central region by a gap in rapidity. 




M (GeV) 

Figure 9: Open squares, AFS data; triangles, tttt cross sections derived from phase shifts of Ochs 
[16]; curve, the elastic cross section with the numerator N{s) replaced by a constant. 

AMP draw attention to structure in the vrvr mass spectum similar to the dip in the tttt S-wave 
elastic cross section. Fig. 9 shows the AFS data as open squares. There is a peak at low tttt 
mass, close to that of the a pole in BES data, but not quite identical, for reasons discussed 
shortly. 

Triangles on Fig. 9 show tttt elastic cross sections derived from Cern-Munich phase shifts; 
this is done by dividing the tttt amplitude of elastic scattering by N{s), leaving only the term 
1/D{s). The curve shows my fit to elastic data after dividing the tttt amplitude by N{el). 
The agreement between the curve and triangles demonstrates that the parametrisation of the 
a reproduces Cern-Munich phase shifts. Results are similar to the AFS data, but there is a 
distinct difference in the vicinity of /o(980). What is clearly evident is constructive interference 
between /q and a immediately below 989 MeV, where EU predicts a zero. Even without fitting, 
one can see that EU will fail. 

A detail is that, up to the KK threshold, one piece of information from each Cern-Munich 
moment is sufficient to determine phase shifts for both S- and P-waves. Above the KK threshold, 
however, the separation of inelasticity parameters 1] and phase shifts 6 cannot be made without 
further assumptions. Just above the KK threshold, the solution of Ochs |1^ has r] =0.6-0.7, 
whereas the BES line-shape for /o(980) demands rj parameters dropping to ~ 0.2 at 1.01 GeV. In 
view of this large discrepancy, predictions from Cern-Munich phase shifts are not shown above 
the KK threshold. 

Since the work of Morgan and Pennington [29], information on both a and /o(980) has 
improved enormously. The a amplitude is taken (within errors) from Ref. [13], using fit (iii) 
given there. The resulting cross section is shown by the chain curve of Fig. 10(a). Although 
this fit follows the general features of the data, it is not accurate for low vrvr masses. Varying 
parameters of /o(980) within errors quoted by BES has negligible effect on the quality of the fit. 
Some additional flexibility is clearly needed in the numerator X{s) of the a amplitude. 

There are two obvious sources. Firstly, AMP point out that central production goes via 
two intermediate Pomerons: PP vrvr, where vr exchange is allowed. This alone does not 
achieve a good flt. The second explanation arises from Regge factors in the production process. 
AMP argue for an s-dependent factor (m^ -|- s) in the numerator arising from the leading Regge 
trajectory. There may be contributions also from a daughter trajectory. 
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Figure 10: Fits to AFS data: (a) full curve isobar model, chain curve using the a amplitude of 
Ref. [15], dashed curve the scaling factor for the final fit. (b) enlargement of the /o(980) region: 
full curve isobar model, chain curve EU fit, dotted EU + freely fitted 1 = 2 S-wave. 



The full curve of Fig. 10(a) shows an isobar model fit using a complex constant for /o(980), 
but without any /o(1370). It uses a scaling factor for the cross section shown by the dashed 
curve. This takes the form of a Gaussian dip: N{prodn) = 1 — Aexp —{s — sq)^. This scaling 
factor gives rise to a slow modulation of the a amplitude over hundreds of MeV. It has only a 
small effect on what is fitted to /o(980). A detail is that a Gaussian mass resolution of 10 MeV, 
quoted by the AFS collaboration [lOj, is folded into each point of the fit; it is significant only 
near the KK threshold. The fit to vrvr data has a of 35.2 for 34 degrees of freedom. 

Consider now the /o(980) mass range. The isobar model fit of Fig. 10(a) requires an /q 
magnitude only ~ 60% of the EU prediction. This fit requires constructive interference between 
a and /o below the /o mass and destructive interference above. It implies the /o phase is (57±7)° 
below the EU prediction. This immediately throws doubt on EU. 

A fit using EU is conspicuously bad. As an example, if the first 29 points up to 988 MeV 
are fitted alone, = 141.7 if fitted by a and /o(980) only. EU predicts an amplitude which 
is almost zero at 989 MeV, whereas the data at 988 MeV are far above this. If the EU fit is 
extended to 1300 MeV including the /o(1370) with the elasticity fitted in [5], is > 10^, because 
the /o(1370) contribution is far too large. Even if the /o(1370) is fitted freely in magnitude, 

= 137 for 37 points. The fit, shown by the chain curve of Fig. 10(b) is particularly bad close 
to the KK threshold, where the dip of elastic scattering is predicted by EU. The fit to KK data 
is also poor, with a of 18 for 5 points. 

The dotted curve of Fig. 10(b) shows the effect of fitting freely an additional contribution 
from the 1 = 2 S-waves: = 59.6. This fails to cure the poor fit. It makes almost no 
difference whether the 1 = 2 amplitude is divided by a factor (s — sa) like the a amplitude. 
The basic difficulty is that the slowly varying 1 = 2 amplitude cannot cure the rapid structure 
due to /o(980). Furthermore, the relative magnitudes of the fitted 1 = 2 amplitude and the EU 
amplitude is 0.46, whereas it is only 0.18 for elastic scattering at 1 GeV. Such a large 1 = 2 
amplitude is implausible. 

Obviously the problem with EU is that the magnitude of the /o amplitude needs to be smaller 
than for elastic scattering. From analyticity, this also requires a difference in phase. That is 
also clear from the absence of the predicted dip at 989 MeV in the data. 

If the phase of the /o(980) amplitude is constrained to the EU value, but relative magnitudes 
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of a and /o are set free, = 144, which is bad. Most of comes from points immediately 
around the KK threshold, showing that the data reject also the phase variation of EU. 



4.1 Proposed Modification to EU 

At this point, one could argue that the mechanism of the production reaction is unknown and 
might generate a phase for /o(980) different to that of the a. This argument is not specific, 
though the isobar model can fit the data well. However, the usual argument for a different 
phase for /q and a is multiple scattering of the pions with spectator particles. In AFS data, 
there is an empty rapidity gap isolating the central region. Remember also that Cern-Munich 
phase shifts are derived in the first instance from data on 7rp ttttp at high momentum and 
small momentum transfer, a similar configuration. 

The conjecture of EU will now be replaced with an alternative ansatz. The treatment of 
production data needs to be able to cope with the case where one resonance amplitude is zero. 
EU does not, since a universal phase equal to elastic scattering requires a production amplitude 
T^^ oc 1/ D'^{s)dI{s). The correct production amplitude should reduce to T'^ when the /o is 
absent. A small /o amplitude should produce a small perturbation to T^. 

Suppose the 2-body tttt amplitude is written as 

Fp, = aT^^ = TTTF^^^ ^ (^TBe''^^% (19) 

where a' and {3 are real. This allows freedom in {3 and includes a phase ^{s) which becomes 
the same as for elastic scattering when /3 — > 1. It is necessary to choose as A the state with the 
larger amplitude on resonance, so that (3 <1. If one could create this tttt system in 'free space', 
the appropriate 2-body unitarity relation below the KK threshold would be 

ImT^^*) - |T(f)|l (20) 

An alternative way of formulating the basic physics (with the same outcome) is in terms of 
the Schwinger-Dyson equation. Instead of the conventional relation 

T^prod ^prod ~l~ ^prod^l^eli (^-^) 

my conjecture is to replace this with 

Tprod ^prod ~l~ ^prod^'-^prod- (^2) 

Here V is the 'potential' generating the final state and G is the propagator. 

For the case of purely elastic scattering, a closed form for ^(s) of (19) may be derived by 
substituting T4 and Tg in the form [e^^^ — l)/2i into (20). After simple cancellations between 
left- and right-hand sides, 

sin(2^ + 6b- 26a) = /3 sin 63, (23) 

or 

2* = 25a -Sb + sm-\p sin 5b). (24) 

If /3 7^ 1, this is a different relation from purely elastic scattering. 

The improvement in the fit is dramatic. Immediately an excellent fit to points below the 
KK threshold is obtained with — 28.4 for 29 points and 24 degrees of freedom. The term 
sm~^ {(3 sin Sb) in (24) differs from Sb by 54°. This is just what is needed to produce the 
interference between a and /o observed in the isobar model fit. 
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Above the KK threshold, it is tempting to satisfy unitarity by introducing a /^-matrix. 
However, the i^-matrix depends on the assumption that the 2-body system is confined to the 
unitary circle, but that is no longer the case in a 3-body situation. 

The fit may be extended above the inelastic threshold by writing tttt and KK amplitudes as 



a 



T^i 1 + 



9i 



9{ , 



Fkk = a[Tf2 + 7r2^2 + (/3T/2 + eT4)e 
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]/P2 

/3 + 



^92 

9{, 



/Pl 



(25) 



(26) 



Eqs. (25) and (26) expose the exphcit dependence of T12 on the ratio r = \j P2I P\\ the experi- 
mental group divides out the phase space p2 and p\ in the tttt and KK channels. As explained 
in Section 2.1, all T have the numerator of elastic scattering replaced by a constant. 

My T12 is defined so as to contain a factor ^ p\P2 and T22 is defined to contain a factor p^- 
With these definitions, the unitarity relations become 



ImTii 
ImTi2 
ImT22 



|Tiip + |ri2p 

ll-'l2 + -'12-' 22 
. ml +-'21-' 12- 



(27) 
(28) 
(29) 



Eqs. (25) and (26) satisfy these relations by construction, except that \Ef and need to be 
constrained to obey (27)- (29). Above the KK threshold, this is done using freely fitted \E' and 
\Ef' for every data point and introducing into a penalty function which applies (27)- (29) with 
3% errors; in practise these constraints are easily satisfied and discrepancies at the end of the 
fit are below the 1% level; this is well below experimental errors. In fact the KK data are not 
very precise, leaving large fiexibility in \E'' above the KK threshold. In other words, the data 
are easy to fit above the KK threshold, but highly definitive below it. 

A detail is that g2 needs to include form factors both below and above the KK threshold, 
such that it falls quite rapidly on both sides of the threshold; formulae are given in [13] . 
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Figure 11: Fits to AFS data: (a) using the revised form of 2-body unitarity, (b) fit to KK data. 
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The best fit with a and /o(980) alone has = 46.4 for 37 points and 32 degrees of freedom. 
The fit is good up to 1.1 GeV, but is inadequate for vrvr data near 1.3 GeV. This may be cured 
straightforwardly by adding a small /o(1370) component. A technical detail is that the /o(1370) 
amplitude is multiplied by a factor exp(2i\I'") and (24) is iterated; the contribution of /o(1370) 
to Ti2 and T22 is negligible. The resulting fit, shown on Fig. 11(a), has /3 = 0.59 ± 0.06 and 
the /o(1370) amplitude is 0.18 times that of the a amplitude at 1.3 GeV. The for '^'^ data is 
28.7 for 29 degrees of freedom. The Omega collaboration reports a significant contribution from 
/o(1370) to their data on central production of tt+tt" [30]. Their fitted mass and width agree 
closely with the line-shape fitted to /o(1370) in Ref. 

Fig. 11(b) shows the fit to AFS K+K' data. A detail here is that the K+K' data of AFS are 
scaled up by a factor 4/3 to allow for isospin Clebsch-Gordan coefficients in tt+tt" and K^K~ 
systems. The KK data constrain the coefficients of T12 and T22 amplitudes. The lowest AFS 
KK point has small acceptance which may have significant systematic uncertainty |31j . 
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Figure 12: The Argand diagram of the amplitude fitting AFS data; masses are marked in GeV. 

Fig. 12 shows the Argand diagram for the fitted amplitude. This illustrates the form of 
Eq. (20). There is a geometrical relation between the imaginary part of the amplitude and its 
modulus squared, but it is a different relation to EU. The /o(980) amplitude is smaller than that 
of the a and their relative phases are different to EU. The same is true of Fig. 5, the Argand 
diagram fitting BES data; there the /o(980) contribution is very small. 

Let us now return to equations (20) and (24) and review their general features. The last 
term of (24), sin~^(/5 sin 5^) — ^ 5_b as /3 1 and 90°. Furthermore it approaches this 

limit non- linearly as /3 — > 1. The interesting point is that two-body elastic scattering emerges 
as a limiting case of a more general relation. Furthermore it has pathological properties as the 
relative magnitude of A and B crosses 1 (or -1). If B becomes larger than A, it is necessary to 
interchange the roles of A and B. As /3 approaches 1 from below, ^ ^ and one recovers the 
standard result of elastic scattering. Further pathological cases arise if /? — ^ — 1 (impossible in 
the 2-body elastic case). Elastic scattering is in fact a very special case. So EU can fail very 
badly as /3 departs from 1. 

As /3 drops from 1, the phase \l/ measured from the origin of the Argand diagram soon 
changes by only a modest amount over the /q. In this case, the isobar model becomes an 
excellent approximation: the /o has its usual dependence on s through its phase shift b and the 
/o amplitude adds vectorially to the a; the isobar model can then fit a constant phase to both 
quite successfully if the term sin~^(/5sin5) is small. 
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In summary, the modified form of EU suggested in (20) and (24) gives a much better fit than 
EU. The fit of Morgan and Pennington [29] using EU requires an additional third-sheet pole at 
M = 978 — i28 MeV. This additional pole cannot be accomodated by BES (pTcir data, which 
require only a second-sheet pole at 998 — il7 MeV and a broad third-sheet pole at 851 — iAlS 
MeV. These two poles have a natural explanation. If the coupling of /o(980) to KK is allowed to 
decrease gradually to zero, leaving other parameters unchanged, the two poles coalesce towards 
the same pole position M = 968 — i82 MeV; it is the effect of p{KK) which moves the second- 
sheet pole to the KK threshold and moves the third-sheet pole away. There is no pole in this 
mass range from the a amplitude. 

A second remark is that all of (T, K, QjQ (980) and /o(980) may be reproduced as a nonet by the 
model of Rupp and van Beveren, where mesons couple to a quark loop [32j. In this model, no 
additional pole like that of Morgan and Pennington appears in the mass range close to the KK 
threshold. These two results suggest that the additional pole is a consequence of the constraint 
of EU. 

The relation (20) is a new conjecture. Are there forseeable snags? The a and /o may mix, 
and this mixing could be different in elastic scattering and production. This mixing would 
alter the apparent width of /o(980) and could induce an additional phase change relative to 
a. Presently there is no indication of any need for mixing. Such mixing will be absent if a 
and /o(980) have strictly orthogonal wave functions, as is plausible for members of the same 
nonet. On Fig. 2, the a is definitely visible in (pinr data. It would not be surprising if (pirir and 
cuTTTT channels filter out orthogonal combinations of a and /q. From Fig. 2, one can estimate 
the relative intensities of a and /q. It is necessary to allow for the mass resolution, since the 
/o amplitude falls extremely rapidly from its peak at 989 MeV, particularly above the KK 
threshold. Doing this, the intensity of a is 4% of /o at the peak, i.e. 20% in amplitde. This is 
marginally higher than the /o signal fitted to ujirn but within the error, supporting the idea of 
orthogonal amplitudes. 



5 How to fit elastic data above the KK threshold 

Many authors use the ii"-matrix to satisfy unitarity for 2-body scattering, e.g. the coupled 
channels vrvr, KK, rji], rji]' and 4n. The popular approach is to add iiT-matrices of all resonances 
appearing in one partial wave. However, if resonances overlap, as a and /o(980) do, the /T-matrix 
poles occur at masses where combined phase shifts happen to go through 90°, 270°, etc., i.e. at 
~ 750 and 1200 MeV. Firstly, an expansion in terms of these poles is problematical for /o(980) 
unless other factors or high powers of s are included. Secondly, the relation between i^'-matrix 
poles and T-matrix poles is obscure. Any one T-matrix pole is built up from all i^'-matrix poles; 
the converse is obvious. The prescription that S-matrices multiply below the inelastic threshold 
does not appear naturally, but has to be enforced by fitting data. 

An attractive alternative can be constructed following the spirit of Aitchison's approach (for 
2-body scattering). Suppose one combines two resonances according to the prescription 

' l-0.5p,p,{KAKB + KBKAh ^ ^ 

Below the KK threshold, this automatically gives the result that phase shifts add. [Further 
resonances may be combined by iterating this prescription]. A nice feature of (30) is that one 
can write 

K,. = (31) 
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using the same mass M as the usual Breit-Wigner denominator. A second attractive feature of 
(31) is that the amphtude continues naturally through the KK threshold, because of the factor 
PiPj in the denominator. 

My own approach in several papers, [HI [5] has been close to this. All diagonal elements of S- 
matrices are multiplied, as proposed here. Magnitudes of off-diagonal elements of the S'-matrix 
need to be calculated from unitarity relations. For example, for a 3-channel system: 

= (1/2)(1 + |533P - \Su? - \S22\')- (32) 

The phase of these off-diagonal elements has been fitted empirically, whereas (30) would predict 
these phases. This approach successfully fits elastic data, tttt KK and r]r] with one proviso: 
a good fit requires inclusion of mixing between a, /o(1370) and /o(1500) [5j, using the formulae 
of Anisovich, Anisovich and Sarantsev |32j; these formulae are the modern equivalent of the 
Breit-Rabi equation of molecular spectroscopy, generalised to include resonance widths. 

6 Conclusions 

Crystal Barrel data have /o(980) and a components with relative magnitudes seriously different 
to those predicted by EU plus analyticity. Furthermore, they are different in ^5*1 and ^Pi anni- 
hilation. The BES data for J/\E' — > cutt+tt" do not reproduce the deep dip of elastic scattering 
at 989 MeV. AFS data likewise do not contain the same dip at this mass. All these results are 
in conflict with equation (14), which is a direct consequence of EU. This shows unambiguously 
that there must be a major flaw in the hypothesis of EU. 

Experimentalists have dealt with this problem by using complex coupling constants for each 
resonance. However, as emphasised in subsection 2.1, the imaginary part of the coupling constant 
introduces into the numerator of the amplitude an s-dependent phase variation which alters the 
universal phase coming from the denominator [1 — ipK{sy\. This destroys the original idea of 
a universal phase. One might as well fit directly in terms of the T-matrix of each individual 
resonance, along the lines outlined in Section 5. The form of the i^'-matrix suggested there 
would eliminate differences between /T-matrix and T-matrix poles, making interpretation of 
results more direct. 

Not all experimentalists adopt the P- vector approach. Some fit directly in terms of individual 
T-matrices for each resonance, including sequential decays from one resonance to a daughter 
with different complex coupling constants for each decay mode. Ascoli and Wyld [M] and Schult 
and Wyld ^35j consider a multiple scattering series of the type R (12)3 1(23), etc, where 

is a 3-body resonance and brackets indicate resonances in two-body sub-systems; this is a 
unitarity effect of a different form to that considered here. 

In view of the failure of EU in the 4 cases considered here, each new set of data should be 
inspected on its merits. 

Let us examine ways of trying to save EU. Firstly, it is possible that unspecified backgrounds 
can be added to the P- vector so as to side-step the conflict. However, analyticity independently 
limits relative magnitudes of /o(980) and a within 12%. Experimental determinations of a 
and K phases in [13] and [T3] constrain phases within ±5°. The probability that unspecified 
backgrounds can evade Eq. (14) to this accuracy in four sets of data is vanishingly small. 
Furthermore, if such backgrounds are introduced, EU loses any predictive power. 

Secondly, Crystal Barrel data and AFS data cannot be fitted with EU whether or not the 
1 = 2 S-wave amplitude is included. So this is not a satisfactory escape route. 

A third likely possibility, applicable to the first three sets of data, is that pions from sigma 
and/or /o(980) rescatter from the spectator particle, leading to a breakdown of EU. Aitchison 
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himself pointed this out. Today, we known that such graphs have magnitudes typically 25% of 
those of the parent processes before the rescattering. This is sufficient to introduce large phase 
changes in some cases, but not all. 

For all of these three sets of data, the isobar model provides an excellent fit. The production 
amplitude is then written F = aANA{el)/ Da + C(BNB{el)/DB with complex and as- In this 
form, no vestige remains of the constraint that S'-matrices must multiply as in elastic scattering. 
Relative magnitudes of q;^ and as can arise from matrix elements coupling the initial state to 
each resonance. For J/^ — > ujtiti and ^Pi pp — > Stt", the /o component is so small that one 
cannot tell whether the phase alone follows EU or not. However, for ^5*0 pp 37r°, the /o signal 
is large enough to rule out this possibility. 

The fourth point is that one would still expect EU to work for AFS data, but it does not. An 
excellent fit may be obtained by replacing the unitarity relation ImT^P) = t^p^T*; by the new 
relation 

ImT^P) = |T(p)|2. (33) 

This corresponds to the relation 

lmF = FT*^P\ (34) 

rather than the commonly used form ImF = FT*i. 

Equations (20) and (24) have pathological behaviour in the vicinity of the 2-body elastic limit 
P = 1. One can now see the basic problem of EU. It attempts to impose on the 3-body system a 
very special behaviour which is narrowly restricted to 2-body scattering. Away from the special 
case 13 — 1, the isobar model works successfully. 

There are two points about the new unitarity relations. Firstly, it was argued in Section 2.1 
that a universal phase in the denominator of the amplitude also requires, via analyticity, almost 
universal magnitudes; the word 'almost' covers the possibility that there may be slowly varying 
form factors or centrifugal barrier factors in production reactions without corresponding changes 
to D{s). If relative magnitudes of resonances differ by large amounts between 2-body scattering 
and production, their relative phases must also change. 

Secondly, the new unitarity relation (33) succeeds quantitatively in accounting for the ob- 
served relative phase between a and /o(980) in central production. That is a non-trivial result. 
The fit to AFS data then requires only two poles in the vicinity of /o(980), in agreement with 
the EES line-shape (as does the isobar model). EU requires an extra pole for which there is 
no obvious explanation. The form of this new unitarity relation is illustrated by the Argand 
diagrams of Figs. 5 and 12. In both, the /o(980) amplitude is small or fairly small and so is 
sin~^(/5sin 5b) of Eq. (19). As a result, the phase ^ measured from the origin of the Argand 
diagram changes rather little over the /q. This is an extra source of phases appearing in the 
isobar model, and has not been appreciated before. 

However, one then needs to ask whether this new relation can be used universally in the isobar 
model. Does it, for example, correctly describe the relative phases of a and /o in ^5*0 — > 37r° data 
and in J/^' — > ^tttt? The answer is no. For these two reactions, the agreement between data and 
equations (20) and (24) improves substantially over EU. However, there are still discrepancies 
with the new unitarity relation of 20 — 30°, which is still significant. It seems likely that 
rescattering of pions from the spectator introduces some additional phases. 

The new unitarity relation needs to be tested elsewhere. A possible testing ground is in Kloe 
data on — >• 77r°7r'', where both a and /o may contribute, but the decay is electromagnetic; the 
small amplitude from (j) p7r° introduces a perturbation of only 4% in amplitude and only in 
well defined parts of the Dalitz plot. 

The remedy which succeeds well in fitting nearly all data is the isobar model, where both 
magnitudes and phases of resonances are both fitted freely. It needs to be emphasised that 
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experimental groups have adopted the flexibihty needed to fit existing data, so their results 
are essentially sound and are not in question. The hypothesis of EU has mostly been adopted 
by theorists for making predictions. Those predictions now need to be viewed with suspicion. 
Although the scalar form factor is well determined for elastic scattering, it is dangerous to 
assume that this form factor is universal and can predict production processes. 
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